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Abstract 



Q,^ ■ The diverse isotopic and elemental signatures produced in different nucle- 

, osynthetic sites are passed on to successive generations of stars. By tracing these 

chemical signatures back through the stellar populations of the Galaxy, it is pos- 
sible to unravel its nucleosynthetic history and even to study stars which are 
now extinct. This review considers recent applications of "stellar genetics" to 
I examine the earliest episodes of nucleosynthesis in the Universe, in Population 

■ III stars and the Big Bang. 

O I Keywords: stars: abundances, stars: Population II, Galaxy: abundances. Galaxy: 

' formation. Galaxy: halo 

> 

^ : 1 Stellar genetics 

Nucleosynthesis in the first few minutes following the Big Bang produced five isotopes 
in significant quantities: ^H, ^H, ^He, ^He, and ^Li. The last of these accounts for less 
than one part per billion. All other elements were produced later, either in stars or as 
a consequence of stellar evolution. 

Stellar nucleosynthesis theory predicts the yields of a wide range of isotopes, but 
in most astronomical spectra only elemental (rather than isotopic) abundances can 
be measured because the wavelength shift between different isotopes of an element is 
usually well below the intrinsic line width set by thermal and pressure broadening. 
Amongst the few exceptions are and (e.g. Linsky et al. 1995; Buries, Kirkman, 
& Tytler 1999), and ^Li and ''Li (Smith, Lambert & Nissen 1993), whose isotope shifts 
are particularly large because of their small atomic masses. Isotope shifts are usually 
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larger for molecular features, and ^^C and ^^C (e.g. Brown & Wallerstein 1989) and 
certain Mg isotopes (e.g. Shetrone 1996) can be distinguished in CH and MgH bands. 
Although isotope shifts become even smaller for higher atomic masses, differences in 
the hyperfine structure of odd and even isotopes may allow constraints to be placed on 
isotope ratios for some elements as heavy as Ba (Magain 1995) and Eu (Hauge 1972). 
Nevertheless, stellar spectral measurements are generally limited to elemental rather 
than isotopic abundances. 

By measuring the time-evolution of abundances coming from a variety of nucle- 
osynthesis sites — see Figure |l] — such as AGB stars, the cores of supernovae, and 
various radial zones in the envelopes of massive stars, it is possible to discover the roles 
of a wide range of stars over the history of the Galaxy. This is the activity I have 
branded "stellar genetics" : using the isotopic and elemental signatures that are passed 
from one stellar generation to the next to trace the particular nucleosynthesis reactions 
that have occurred in the past, and hence also to trace the sites of those reactions even 
when these may no longer exist. For example, low-metallicity asymptotic giant branch 
(AGB) stars have all now become white dwarfs, but it is nevertheless possible to study 
the internal structure of those AGB stars via the material they have dumped onto the 
surface of longer-lived companions (Ryan et al. 2001a). 
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Figure 1: Isotopic abundances arranged by atomic-mass number A. The legend indi- 
cates the main nucleosynthesis pathway for each species. (Based on data from Lang 
1980, Table 38). 

The supernovae of massive stars are important early contributors to the enrichment 
of the Galaxy. Stars whose main-sequence mass exceeds 8-10 Mq are expected to 
become core collapse (type II) supernovae on timescales < 50 Myr, as short as ~1 Myr 
for stars of a few tens of solar masses. SN II produce of order ten times more oxygen 
than iron (by mass), whereas SN la, with lower-mass progenitors and hence longer time 
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lags, produce of order seven times more iron than oxygen (e.g. Timmes, Woosley & 
Weaver 1995). Consequently, the abundance ratios of oxygen and the a-elementsQ to 
iron are higher in Galactic halo stars than in the Sun (e.g. Tinsley 1979, Norris, Ryan 
& Beers 2001). 

Although supernovae have been widely modelled, the role of stellar winds in remov- 
ing processed material from highly luminous stars has not yet been widely investigated 
for the early enrichment of the halo, presumably because of uncertainties in the treat- 
ment of mass loss and its lower efficiency in metal-poor environments. However, the 
unexplained high levels of carbon enrichment in the early Galaxy, where perhaps 10% 
of stars have carbon-to-iron ratios up to 100 times solar (Norris, Ryan & Beers 1997a; 
Barbuy et al. 1997; Rossi, Beers & Sneden 1999) suggest that more attention should 
be devoted to investigating this means of enriching the early interstellar medium. This 
is not to say that stellar winds are the only way of effecting large carbon excesses; 
indeed, the carbon-rich stars show a range of other signatures including s-process en- 
hancements (Norris et al. 1997a), r-process enhancements (Cowan et al. 1995), and 
normal neutron-capture-element abundances (Norris, Ryan & Beers 1997c, Aoki et al 
2002). Nevertheless, one may speculate that mass-loss may help produce high carbon- 
to-iron ratios without producing neutron-capture elements if the ratio of black-hole 
remnants to neutron star remnants is greater in metal-poor than metal-rich popula- 
tions, as it might be since metal-poor stars are more compact and hence sit deeper 
in their gravitational potential wells. This might manifest itself in a lowering of the 
progenitor mass which corresponds to the division between (low-mass) neutron-star 
remnants and (high-mass) black-hole remnants in metal poor populations. (This un- 
certain value is sometimes set in the range 30-50 Mq in Galactic chemical evolution 
calculations; see also Prantzos (1994).) This would allow stellar winds to enrich the 
interstellar medium in carbon, yet avoid an associated enrichment in iron due to fall 
back onto a neutron star. 

Stars below the mass limit for SN II become AGB stars which are responsible 
for most Galactic nucleosynthesis of carbon and s-process elements, but because of 
their lower masses they enrich the Galaxy after the first SN II. Models of 8-10 
supernovae as the source of r-process nuclei suggest that these appeared as the Galaxy 
reached [Fe/H] ~ —3.0 (e.g. Mathews & Cowan 1990), which implies that AGB stars 
will enrich the halo only once this metallicity has been reached. This makes AGB 
stars unsuitable as the source of the carbon excess in lower metallicity stars. Also, 
the s-process contribution is insignificant until an even higher metallicity is reached, 
because of the need for pre-existing seed nuclei (Spite & Spite 1978; Truran 1981). The 
numerical model of Travaglio et al. (1999) suggests that the s-process appears only at 
[Fe/H] —1, though Pagel & Tautvaisiene (1997) argue for its appearance as low as 
[Fe/H] = -2.5. 

Although the science of stellar spectroscopy is more than a century old, many of the 
oldest stars in the Galactic halo have been discovered only in the last decade, as a result 
of the Beers, Preston & Shectman (1985, 1992) survey and its successors (e.g. Norris, 

^a-elements are those relatively abundant light elements whose nuclei may be viewed as multiples 
of He nuclei, resulting either from He-burning or further fusion of He-burning products. 
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Ryan & Beers 1998; AUende Prieto et al. 2000). Furthermore, detailed spectroscopic 
study of such stars requires echelle spectrographs on 4-10 m-class telescopes equipped 
with efficient electronic detectors. The remainder of this review highlights a number 
of issues concerning the very earliest stages of chemical evolution, namely Big Bang 
nucleosynthesis and the lack of identifiable Population 111 stars. Recent reviews of the 
evolution of other elements in the early phases of the Galaxy can be found in Ryan 
(2001) and Ryan et al. (2001b). 

2 The elusive Population III 

The term "Population III" means different things to different people. For some it 
refers to very metal-poor stars, whereupon the distinction between Pop 111 and Pop 11 
becomes blurred. More usefully it can be applied to the first generation of stars that 
formed from material enriched by Big Bang nucleosynthesis alone, prior to stellar 
nucleosynthesis occurring. This is the definition that will be adopted here. The massive 
stars of this ffist stellar generation were probably responsible for the re-ionization of 
the universe, possibly at redshifts around 2; ~ 10 (Tumlinson & Shull 2000), before the 
formation of Galactic halo stars at a redshift 2; ~ 5 (e.g. Edmunds & Phillipps 1997). 

It might be supposed that searches for extremely metal-poor stars should encounter 
some having Population III composition, but none has been recognised. In fact, few 
would ever have existed. Recall that a single 25 M© supernova is sufficient to en- 
rich a 10® M0 primordial gas cloud to [Fe/H] ~ —3.5 (Ryan, Norris & Besscll 1991), 
demonstrating that very little stellar evolution and nucleosynthesis is required to sig- 
nificantly enrich the interstellar medium. Furthermore, such supernovae will appear on 
timescales of only a few million years. Pop III stars will be a very minor constituent 
of the Galaxy. 

A more instructive question to ask is whether there is a lower limit on the metallicity 
of Pop 11 stars. The answer appears to be "yes". Despite the systematic surveys cited 
above, no star more metal- deficient than CD— 38° 245 (Bessell & Norris 1984) has been 
discovered, though a few of comparable metallicity have been found. Although a simple 
closed-box model for galactic chemical evolution is a poor model for the Galactic disk, it 
docs a surprisingly good job of fitting the Galactic halo (Hartwick 1976; Ryan & Norris 
1991a) and predicts a factor of ten fewer stars for each factor of ten lower metallicity. 
However, the good fit breaks down once [Fe/H] falls to —4, and a substantial deficit 
of stars becomes apparent. Based on current halo samples, the simple model suggests 
that 11 stars with [Fe/H] < —4 should have been detected, whereas at most two are 
known (Norris 1999). The Poisson probability of this occurring by chance is just 0.12%, 
strongly suggesting the deficit is real. 

A number of potential explanations for the absence of the more metal poor stars 
can be proposed: 

• The lowest metallicity (Pop III?) stars probably formed before the Galaxy had col- 
lapsed to the size of the halo, and hence the missing stars may have a greater radial 
scale than the halo. 

• The initial mass function (IMF) may have had a lower limit above 0.8 Mq, so no 
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main-sequence or red giant branch Pop III stars remain after ~15 Gyr (e.g. Nakamura 
& Umemura 1999). 

• Since the pre-main-sequence contraction time of 0.8 stars is ~ 10^ yr whereas 
the first supernovae would have exploded after only ~ 10^ yr, some contamination of 
low-mass Pop III stars may have occurred. However, since low-mass stars will have al- 
ready contracted to the base of their Hayashi tracks and undergone deuterium burning 
on a ~ 10^ yr timescale, it is less likely that contamination from a supernova would be 
very effective. 

• Pop III stars may have accreted sufficient material on passages through the Galactic 
disk that their thin outer convcctive zones are now contaminated by heavier elements 
(e.g. Yoshii 1981). However, such accretion could not hide Pop HI giants, whose 
greater convection would give them lower metallicities than dwarfs (which is not ob- 
served) , and their heavy element abundance ratios would resemble those of the accreted 
material, but halo abundance ratios are different from that of disk material (Ryan & 
Norris 1991b). 

So, Pop HI stars are undetected, but that does not mean we cannot study them. 
Thanks to stellar genetics, the yields of the first Pop III supernovae are observable as 
the post-primordial composition of Pop II stars. Comparisons of the abundances of 
extremely-metal-poor halo stars with the yields predicted for Pop HI supernovae (e.g. 
Woosley & Weaver 1995) search for concordance in the abundance ratios, with the aim 
of revealing the mass, supernova-energy ranges, and possibly eventually the rotation 
of Pop III stars, thus identifying which types of stars have contributed to the early 
chemical enrichment of the Galaxy (e.g. Norris et al. 2001). 

3 Primordial isotopes 
3.1 ^He 

Most Pop II stars are not hot enough to excite He spectral hnes, and those that are 
show the effects of diffusion which, in the absence of convective mixing, causes elements 
to stratify in the stellar atmosphere. To find the primordial ^He composition, it is 
necessary to measure H II regions in galaxies that are less chemically enriched than 
the Milky Way. By measuring the strengths of emission lines of hydrogen, helium, 
and various metals (especially oxygen and nitrogen), it is possible to derive relative 
abundances (e.g. Pagel et al. 1992). 

The evolution of the helium mass fraction Y with metallicity mass fraction Z is 
revealed in a graph of Y vs Z, which gives a straight line graph of the form Y — 
Yp + {dY/dZ)Z. The primordial helium abundance Yp is simply the intercept. 

Challenges in deriving the primordial helium mass fraction Yp include: 

• Helium production does not always accompany metal production at the same rate. 
That is, dY/dZ is not the same for all objects. Wolf-Rayct stars enrich H II regions 
more in He than the metals suggest, and bias Yp upwards. Fortunately, Wolf-Rayet 
stars have a spectral signature which allows the most troublesome cases to be detected 
(e.g Pagel et al. 1992), but some bias may remain at a low level. 
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• H II regions are visible only because the gas is excited by the UV thermal radiation 
from hot stars embedded in the clouds. The absorption lines in the stellar spectra lit- 
erally undermine the emission line intensities, resulting in artificially weak He emission 
lines, from which low values of Yp would be derived (e.g. Olive & Skillman 2001). 

• It is often assumed that only UV thermal radiation from the hot stars ionises and 
excites the nebula gas, but there may also be some excitation by collisions between gas 
atoms or by re-absorption of the emitted radiation. The effect of ionisation correction 
factors is unclear (Ballantyne, Ferland & Martin 2000; Viegas, Gruenwald, & Steigman 
2000). 

An illustration of how these factors have influenced primordial He estimates is given 
in Figure |^. All values are close to Yp = 0.24, but since the mid-1990s two streams 
of values have been obtained: a "low" value near 0.230, and a "high" value near 
0.245. Moreover, the stated error bars are considerably smaller than this difference! 
The impact of these different values on the derived baryon density of the universe is 
discussed later. 
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Figure 2: Measurements of the primordial helium abundance, 1979-1998. Horizontal 
points and bars give the inferred values, with references in the right-hand column. 



3.2 ^He 

^He is presently almost unusable as a diagnostic of Big Bang nucleosynthesis. In the 
solar system, direct measurements of ^He atoms in the solar wind and meteorites can 
be made, but all of this material has passed through numerous stellar generations. The 
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only other diagnostic is the 3.4 cm hyperfine transition of the ^He II ion, analogous 
to the 21 cm hyperfine transition of I and the 92 cm hyperfine transition of I, 
which can be observed in H II regions and planetary nebulae. 

To derive an abundance from the observations, an accurate model of the H II region 
or planetary nebula is needed. There is reason to doubt the suitability of current 
models, because different H II regions show a range of ^He abundances that correlate 
with the mass of the H II region, suggesting that the ^He has been depleted locally 
(Rood et al. 1995). 

There are also large uncertainties in the interpretation of ^He abundances. Stars are 
expected to both produce and destroy ^He at different stages of their life. Until very 
recently, the expectation was that 1-2 Mq stars would be net producers of ^He while 
more massive stars would destroy it, but that the overall effect would be an increase 
in ^He with metallicity. This was a problem, because it predicted much more ^He than 
is observed. This may be resolved by recent developments in stellar evolution theory, 
through the proposal by Boothroyd & Sackmann (1999) and Sackmann & Boothroyd 
(1999) that mixing of material occurs deeper in stars than would be expected from 
standard convective theory, and that this leads to destruction of ^He even in the lower 
mass stars previously believed to be net ^He producers, making stars overall destroyers 
of this primordial isotope. 

Given these significant problems, ^He currently does not constrain primordial nucle- 
osynthesis strongly. Bania, Rood & Balser (2001) describe their measurement ^He/H 
= (1.79±0.65)xl0"^ as a "reasonable approximation of the primordial value," while 
Bell (2000) reached the contrasting conclusion that "the upper hmit for the relative 
abundance of primordial ^He ... is ^He''"/H+ < 2 x 10~^." 

3.3 2h 

Deuterium is destroyed in stars, burning as d + p ^^He + 7. As there is no net 
source of deuterium in stars, its abundance has decreased steadily since the Big Bang, 
and any value measured today is a lower limit on the primordial value. Accurate 
measurements of the deuterium abundance can be made in solar system bodies and 
the interstellar medium of our Galaxy. The interstellar "Da" line, which is displaced 
from the "Lyman-a" line of at 1216 A by —0.33 A, gives D/H = 1.6x10"^, which 
is therefore a firm lower limit on the primordial value (Linsky et al. 1995). Models of 
the chemical evolution of the Galaxy can be applied to estimate how much deuterium 
has been destroyed over its lifetime, but this adds more uncertainty. 

Much effort since the mid-1990's has gone into trying to measure deuterium in 
the interstellar medium of other galaxies, particularly in damped Lyman-a (DLA) 
absorption systems superimposed on quasar spectra at high redshift. The look-back 
times of objects at high redshifts offer the opportunity of looking directly into the 
past, and making measurements of the deuterium content "back then". Such values 
are expected to be closer to the primordial one than values measured in the interstellar 
medium of our own Galaxy. 

Because the abundance of deuterium is so much lower than that of hydrogen, large 
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column densities of hydrogen are required before deuterium can be detected. Damped 
Lyman-a clouds have high column densities and are believed to be early galactic gas. 
Debate continues as to their nature: are they collapsed into flattened disks as in a 
mature galaxy like the Milky Way, or still distended like the halo of the Galaxy? The 
restriction to using high column-density clouds is the source of several difficulties in 
obtaining reliable D/H measurements at high redshift. 

• Firstly, there is confusion over whether any single absorption line is due to deuterium 
or to another hydrogen line at a slightly lower redshift (e.g. Buries et al. 1999). 

• Secondly, although it may be possible to measure the deuterium column density quite 
well, we also need to know the hydrogen column density if we are to obtain the ratio 
D/H. However, since the Lyman-a line is saturated (has absorbed all of the light), it 
is very difficult to measure its column density accurately. The H column density has 
to be obtained from other features in the spectrum. 

Figure ^ shows measurements towards a range of quasars over 1994-1999. A con- 
siderable range of values has been obtained, including "low" values around D/H = 
(2-5)xl0-^ and "high" values around D/N = (15-30) xlO'^ 
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Figure 3: Measurements of the deuterium abundance in quasar absorption line systems, 
1994-1999. (See Figure |^ for the layout.) The value measured in the interstellar 
medium (ISM) of the Milky Way is shown at the bottom for comparison. Below the 
figure are recorded the quasar names, the redshift of the measured absorption system 
(not the redshift of the quasar, Zem), and the metallicity Z on a logarithmic scale 
relative to the Sun, where [Z/H] = logio (Z/Zsun)- 
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The high and low values have very different implications. As any real measurement 
of deuterium sets a lower limit on the primordial abundance, if any of the "high" 
measurements is correct, then the primordial value must be (at least!) an order of 
magnitude higher than that inferred from the Galactic ISM. However, the difficulty 
of knowing whether absorption just shortward of the damped Lyman-a line is really 
due to deuterium weakens the arguments for very high D/H ratios. "Low" values 
have been measured in several systems, in particular the ones with low metallicity 
[Z/H] ~ —3.0, —2.9, and —2.8. For the "high" values to be correct, deuterium would 
have to be destroyed very quickly in the low metallicity systems despite their low 
metallicities indicating very little stellar processing. This seems unlikely (Tosi et al. 
1998; Prantzos & Ishimaru 2001). It is more likely that the lowest metallicity systems 
are more reliable measurements, and that the "high" values are misleading due to the 
presence of additional absorption lines. Data are still being taken; it is perhaps a 
matter of opinion whether the issue has been settled, but it is probable that the "low" 
deuterium measurements will win out over the "high" ones once the dust settles. (See 
references to Figure H.) 

3.4 ^Li 

Lithium is measurable in the spectra of main-sequence halo dwarfs. ^Li is much less 
fragile than ^H, and as the outer layers of solar-mass metal-poor dwarfs do not mix 
to great depths, the surface Li in these stars is widely believed to be little changed 
since they formed ~13 Gyr ago. The main evidence in support of this is empirical: 
most solar-temperature, metal-poor dwarfs have almost the same Li abundance despite 
having different amounts of other elements and different masses (e.g. Spite & Spite 
1982). It is unlikely that all such stars could deplete or supplement their Big Bang Li 
complement and still end up with the same value after ~14 Gyrs. Theoretical models 
of stars differ in the amount of Li they preserve, and it is just conceivable that some 
mechanisms of Li destruction might conspire to deplete all stars by the same amount. 
However, the observations can tolerate little Li depletion even over the long lifetimes 
of these stars. For processing by rotationally-induced mixing, Ryan, Norris, & Beers 
(1999) and Pinsonneault et al. (2001) infer limits < 0.1 and ~0.2 dex respectively. 
The difference in these values depends on (A) the frequency of the few stars with Li 
abundances significantly lower than the rest, (B) the size distribution of their defi- 
ciencies, and (C) whether the deviations are in fact due to the proposed mechanism 
(rotationally-induced mixing) and not some other. Small amounts of Li production 
may also be present in the stars that formed more recently from nucleosynthetically 
enriched material. This can be tracked analogously to the way He production is tracked 
by dy/dZ. The main uncertainty remaining for Li is how accurately we know the tem- 
peratures of the stellar atmospheres, and therefore how well we can compute the Li 
abundances from their spectra. The estimate of the primordial Li abundance by Ryan 
et al. (2000) is n(Li)/n(H) = 1.23l|]:| x IQ-^o. 

Note, however, that even though the primordial Li abundance may be reasonably 
well known, there are several significant deficiencies in our understanding of Li pro- 
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cessing in Pop II and Pop I stars. About 7% of Pop II stars which would be expected 
to have normal Li abundances show extensive (essentially complete) destruction of the 
element. Investigations of other elements in these stars (Norris et al. 1997b; Ryan, 
Norris & Beers 1998) show heterogeneous behaviour and indicate that no single char- 
acteristic besides Li deficiency unites these stars. Ryan et al. (2001c) have suggested 
they may be blue-stragglers-to-be, sub-turnoff-mass products of classical blue-straggler 
formation mechanisms. There is also an unexplained contrast between the highly uni- 
form Li abundances measured in halo field stars (Ryan et al. 1999) and the large 
spread of measurements for subgiants in the halo globular cluster M92 (Boesgaard et 
al. 1998). 

3.5 Putting it all together 

We now put the observational constraints together and examine the concordance with 
theory. The standard Big Bang model assumes a uniform density and predicts the 
yields of the primordial isotopes for different values of that density.0 The predicted 
yield of ^He is close to 24% for a wide range of densities, and the fact that this is 
so close to the real universe is a remarkable accomplishment of the model. It is even 
more remarkable that for and ^Li the abundances are also in accord. For three 
isotopes spanning a range of 1 billion in abundance, the model is broadly in accord 
with observations. 

Figure ^ presents the situation in more detail, and illustrates the different values of 
the derived ^He and abundances, and the ^Li value. No ^He observations are shown. 
The predicted yield of Li is not monotonic with density, as a second nucleosynthesis 
pathway becomes effective for rj > 2x 10^^°. The well-defined Li abundance is therefore 
consistent with a wide range of density values. Two observational results are given for 
each of ^He and ^H, the "high" and "low" values discussed above. Corrections for 
systematic errors in the ^He analysis favour "high" values, and the study of more 
damped Lyman-a systems favour the "lower" range for D, giving a value for rj ^ 
5 X lO^^*'. (The slightly lower 95% confidence interval given by Ryan et al. (2000), 
7] = (1.7 — 3.9) X 10~^°, was based on a primordial ^He abundance intermediate between 
the two observational values shown here.) 

An independent measurement of the universal baryon density is possible from the 
power spectrum of the cosmic microwave background radiation (CMBR). The higher- 
order fluctuations above the main (~ 1°) peak are influenced by the propagation of 
sound waves, which reflects the baryon density. Measurements published in 2000 giving 
?7 = 10 X 10~^° raised the spectre of incompatible Big Bang nucleosynthesis (BBN) and 
CMBR measurements of the baryon density (Hogan 2000). Fortunately, the earlier 
values have been superseded with the addition of more data such that the CMBR 
measurements now suggest a 95% confidence interval r] = (3.4 — 7.9) x lO"^*', which is 
certainly compatible with the BBN values cited above (Netterfield et al. 2001). 

^The density may be expressed variously as the baryon-to-photon ratio, ry, or the baryon density 
as a fraction of the critical density, ^b, where ^ 0.014 x x ( 50 km s-^ Mpc-i -^2_ 
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For an assumed Hubble constant of 63 km s~^ Mpc~^, the value rj = 5 x 10^^° 
corresponds to Qb = 0.043. Not only is this well below the closure density, it is also 
a factor of ten below the matter density Qm — 0.35 inferred from other cosmological 
studies (Efstathiou et al. 1999), pointing to the continuing need for non-baryonic dark 
matter. 
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Figure 4: Solid curves: Calculations of yields of Big Bang nucleosynthesis for ^He, ^H, 
the sum of and ^He, and ^Li, for different assumed baryon-to-photon ratios. Long 
dash curves: 95% confidence range for the nucleosynthesis calculations, primarily re- 
flecting uncertainties in the nuclear reaction rates. Data symbols: current observational 
values, split for ^He and into "high" and "low" values to reflect current systematic 
uncertainties. 
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